A (2+1) dimensional integrable spin model: 
Geometrical and gauge equivalent counterpart, solitons 
and localized coherent structures 



R. MYRZAKULOV"'*, S. VIJAYALAKSHMI 6 , 
G. N. NUGMANOVA" and M. LAKSHMANAN 6 't 
l High Energy Physics Institute, National Academy of Sciences, 
480082, Alma-Ata-82, Kazakstan 
^Centre for Nonlinear Dynamics, Department of Physics, 
Bharathidasan University, Tiruchirapalli 620 024, 

India 



Abstract 

A non-isospectral (2+1) dimensional integrable spin equation is investigated. It is 
shown that its geometrical and gauge equivalent counterparts is the (2+1) dimensional 
nonlinear Schrodinger equation introduced by Zakharov and studied recently by Stra- 
chan. Using a Hirota bilinearised form, line and curved soliton solutions are obtained. 
Using certain freedom (arbitrariness) in the solutions of the bilinearised equation, ex- 
ponentially localized dromion-like solutions for the potential is found. Also, breaking 
soliton solutions (for the spin variables) of the shock wave type and algebraically lo- 
calized nature are constructed. 
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The nature of (1+1) dimensional integrable systems is now well understood [1]. On the 
other hand, examples of (2+1) dimensional integrable equations solvable by Inverse Scat- 
tering Transform method (1ST) are fewer in number and such systems are being actively 
investigated from different points of view at present [1,2]. An interesting subclass of inte- 
grable systems, useful both from the mathematical and physical points of view, is the set 
of integrable spin systems. Since the identification of the first integrable spin model twenty 
years ago, namely the continuum isotropic Heisenberg spin systems [3,4], several other in- 
tegrable spin systems in (1+1) dimensions have been identified and investigated (see for 
example, refs.[5,6]) through geometrical and gauge equivalence concepts and 1ST method. 

Again in (2+1) dimensions, only a small number of integrable spin systems are known, 
among which the Ishimori equation (IE) is the most prominent one admitting different 
kinds of spin excitations such as solitons, vortices, dromions and so on [2]. Since the spin 
equations can also be considered as the nonrelativistic version of 0(3) sigma models, which 
in turn have wide physical ramifications in (2+1) dimensions in problems such as high 
T c superconductors and quark confinement, construction of any integrable spin model in 
higher dimensions assumes considerable significance. Moreover, since the spin systems have 
a natural connection with nonlinear Schrodinger family of equations in general, it is relevent 
from the soliton theory point of view to establish such connections in (2+1) dimensional 
cases also (as in the case of IE and Davey -Stewartson equation). 

Recently a new family of integrable and nonintegrable (but admitting exact solitary wave 
solutions) (N+l) dimensional (N = 1,2) classical spin models was proposed in refs[7-9]. One 
interesting nonlocal integrable spin model is (the so called M-I equation) [7] 



S t = {SAS y + uS} 





(16) 



Here, subscripts stand for partial derivatives, S = (Si, S 2 , S3) and S 2 = Sf + r 2 (Sf + Sf) = 1, 
r 2 = ±1. As in the case of IE, the quantity 
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q = dxdyS ■ (S x A S y ) 




(2) 



may be called the topological charge. Both eq.(l) as well as IE in the (1+1) dimensional case 
reduce to one and the same equation - the well known (1+1) dimensional isotropic classical 
continuous Heisenberg ferromagnet model [3]. The Lax representation of eq.(l) was given in 
ref. [7] and some of its properties were studied in refs.[8-ll]. The aim of this letter is to find 
the equivalent (geometrical and gauge) counterpart namely the (2+1) dimensional nonlinear 
Schrodinger equation(NLSE) and to obtain physically interesting solutions such as line and 
curved solitons, exponentially localized and breaking solitons of eq.(l). 

To begin with, let us find the geometrically equivalent counterpart of eq.(l), for r 2 = 1, 
that is, when S 2 = S 2 + S 2 + S| = 1. For this purpose, we will extend the geometrical 
method applicable to (1+1) dimensional systems suitably to the (2+1) dimensional case. 
We associate a moving space curve parametrised by the arclength x, and endowed with 
an additional coordinate y with the spin system[3,15,16]. Then the Serret-Frenet equation 
associated with the curve has the form 



and e*j's, i = 1,2,3 form the orthogonal trihedral. Mapping the spin on the unit tangent 
vector 



D Ae* 



(3a) 



where 



D = re*i + /te*3 



(36) 



S(x,y,t) = ei, 



(4) 



the curvature and the torsion are given by 



K(x,y,t) = 



(5a) 



r(x,y,t) = k 2 S- (S x A S xx ). 



(56) 
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Due to the orthonormality nature of the trihedral, e it .ei = 0, e iy .ei = 0, i,j = 1,2,3 and 
using the compatibility condition e ixy = ei yx , we find the equation for the y-part 



where 7 = (71, 72, 73) and 



(6) 



7! = « + c\ Vy, 



(7a) 



72 = — 



u 3 



K 



(76) 



73 = 9- fa - ^) . (7c) 

Now, from eq.(l) and using eqs.(3) and (6), we can easily find the time evolution of the 
trihedral 



e it = oj A e i: 



(8) 



with 



Q = (wx, w 2 , w 3 ) = ( ^ - rd s \, -Ky, -nd x 1 T y ) . 



(9) 



Ultimately, the compatibility condition e ixt = e^, which is also consistent with the 
relation e iyt = e ity , i = 1, 2, 3 yields the following evolution equations for the curvature and 
torsion 



Kt — {^KT) y K,xd x Ty, 



(10a) 



n 



^1 _ r g-i r 

IU X ly 



H - KjKjy. 



(106) 



On making the complex transformation [3], 



K{x,y,t) 
V{x,y,t) = exp 



/X 
t(x', y, t)dx' 
-00 



(11) 
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the set of equations (10) reduce to the following (2+1) dimensional nonlinear Schrodinger 
equation(NLSE) 



= i>x y + r*Vi/>, (12a) 

V x = 2d y \ i> | 2 . (126) 

Here, r 2 = +1, that is, we have the attractive type NLSE (The case r 2 = —1 corresponds 
to the repulsive case). Eq.(12) was originally introduced by Zakharov[12] and was recently 
rederived by Strachan (for r 2 = +1)[13]. Its Painleve property and some exact solutions 
were also obtained [14]. N- soliton solutions of eq.(12) for both the cases (r 2 = ±1) can be 
found in ref. [17]. Here, we have proved that eq.(12) is equivalent to eq.(l) in the geometrical 
sense. 

Next, it is always of interest to note that eqs. (1) and (12) are also gauge equivalent in 
the sense of Zakharov and Takhtajan[18]. To this end, we write the Lax representation of 
eq.(l)[7], 

<f>ix = Ui<f)i, (13a) 



<t>it = V 1 (f) 1 + \<f> ly , (136) 



where 



V 1 = ^([S,S y ] + 2vuS), S ± = S 1 ±iS 2 . (146) 

Here, A is the eigen value parameter which satisfies the following equations of Riemann wave 
type 

A t = XXy (15) 
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It means that for solving eq.(l), we must in general use the non-isospectral 1ST. To obtain 
gauge equivalent counterpart of eq.(l), in the usual way we consider the following gauge 
transformation 

0i = <rV 2 , (i6) 

where g(x,y,t) and (f>2(x,y,t, A) are temporarily arbitrary matrix functions. Substituting 
eq.(16) into eq.(13), after some algebra we get the following system of linear equations for 

02 

<p2x = ^202, (17a) 



4>2t = V 2 <f>2 + A0 



2y 



(176) 



with 



2 \ -r 2 (b 



(18a) 



(VI 

V 2 = icr 3 { — -Gy), I = diag(l,l), 



(186) 



V = 2d-% (| i> | 2 ) (18c) 

The compatibility condition of eq.(17) with (15) becomes (12), that is, eq.(l) and eq.(12) 
are gauge equivalent to each other. The above transformation is in fact reversible and we 
can similarly prove that eq.(12) is gauge equivalent to eq.(l). 

The integrable eq.(l) allows an infinite number of integrals of motion. The first two 
conservation laws are 



(3), + ^ 



g2^-i i s ■ S x A S x 



s?. 



= 0, 



+ 



(19a) 
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S ■ (S x A S xx ^j 



+ -d x 
t 2 



(S x ) X (S%)y 



S-(S x AS xx )d^ (S-(S x AS xx )) y 



+d y { 



(^+ 2(g ' ( i Ag " ))2 -2(a -4g 



C2 



> = o 



(196) 



and so on. Next, we present some important formulae which are just consequences of geo- 
metrical/gauge equivalence of eqs. (1) and (12). We have 



tr(S 2 x ) = 8 | V ?= 2S : 



(20a) 



In a similar manner we find that 



-2iS ■ (S x A S xx ) = tr{S x SS xx ) = 4(^^ - (206) 

These relations are obviously equivalent to eq.(ll). One may note that these are of the same 
form as in the case of (1 + 1) dimensional Heisenberg chain [3]. 

Now, we wish to find a class of exact solutions of eq.(l), such as line and curved solitons 
as well as exponentially localized solutions. Introducing the stereographic variable 



S+ = S 1 + iS 2 



2oo 



1+ I U) 



1- I u 
1+ I to 



(21) 



eq.(l) takes the form 



• / \ 1uJ*LO x U3y 

I (U t ~ WjJ x ) + U xy ~ T~ I9N = °' 



(1+M 2 ) 



2i(u x u* - u* x u) y ) 

Ux+ (l+~FF? = 



(22a) 



(226) 



On writing 



g_ 

/' 



(23) 



where g and / are complex valued functions, and after using the Hirota's D-operators, eq.(22) 
becomes 
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(iD t -D x D y )(f*og) = 0, 



(24a) 



(iD t -D x D y )(f*of-g*og) = 0, 



(246) 



D x (f*of + g*og)=0, (24c) 

while the potential u(x,y,t) is 

uix.y.t) — —i — ^ . (24a) 

V ' f*of + g*og V ; 

— * 

In terms of g and /, the spin field S takes the form 



Q+ — 2 / j q _ I / I I 9 I /r,rN 

D ~~ I -f 12 i I 12' ' 53 ~~ I -f 12 i I 12 ^°' ) 

I / 1 2 + I g r I / 1 + \9 1 

and for u(x, y, t) we get the following formula (using the properties of .D-operators) 

u x (x,y,t) = -2id 2 xy \n(\f\ 2 + \g\ 2 ). (26) 

The construction of the solutions to the M-I equation (1) now becomes standard. One 
expands the functions / and g as a power series in the arbitrary parameter e, 



n=0 n=l 

Substituting these expansions into (24 a,b,c) and equating the coefficients of e n yields 



[idt + d x d y ] g 2n +i = - J2 D '(fL ° 92m+i), (28a) 

k+m=n 



id t - d 2 xy 



(fL - /an) = D' ( J2 9*2n 1+ l°92n 2+ l- £ / 2 * mi ° / 2m2 ] , (286) 

ni+ri2=n— 1 mi+m,2=ri 



^(/ 2 *n - /2n) = D x [ ]T 9*2n 1+ l° 92n 2+ l~ ^ fL^f**}, (28c) 
ini+n2=n— 1 ni+n,2=n 
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with D' = iD t — D x D y , f = 0. In order to construct exact N-soliton (line and curved) 
solutions (N-SS) of eq.(l), we make the ansatz 

N 

g 1 = J2 ex PXj, Xj = X j x + m j (y,t) + c j . (29) 

We note here the important fact that m,j (y, t) is an arbitrary complex function of (y, t) of 
the form (see eq.(28a)) 

mj(y,t) = rrijip), p = y + i\jt, (30) 

where Xj is an arbitrary complex parameter. As an example, we write the forms of g and / 
for N=l as 

£i = expxi, f 2 = exp2(xii? + 4>), (31) 

where 

Xi = Xir + i\u, Xi = rj + i£, m 1 = m 1R (p) + im u (p), Xir = V% + m 1R (p) + c 1R , 

—A 2 

Xu = & + mu(p) + cu, c = ln(2rj/ A*), exp 2ip = — — — 

(Ax + AJ 2 

m lR (p) = Remi(p),mu(p) = lmm x {p). (32) 
The corresponding 1-SS of eq.(l) takes the form 
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Ss(x,y,t) = l ^—jsech 2 XiR, (33a) 



2r] 

while for the potential 



S + (x, y, t) = ^ 2 [if - rjtanhxi R ] sechxm, (336) 



y, t) = (f "lifl - Wu) sec/l Xir (33c) 

_|_ ^ 

in which the prime has been used to denote the differentiation with respect to the real part 
of the arguement. We note that the 1-SS(33) depends on two arbitrary functions mi R (p) 
and mu(p) as in the case of some other (2+1) dimensional integrable equations. 
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Now for the particular choice, 



mi = ki(y + i\it), 

where k\ is a complex constant, in eq.(33), we get the usual line 1-SS of eq.(l). For fixed 

— * 

(y,t), it follows from (33) that S — > (0,0, 1) as x — > ±oo and the wavefront itself is defined 
by the equation xir —T]x + m\R{p) + cm = r/x + k\y — + ci# = 0. For other choices of mi, 
we can obtain more general solutions. Particularly, we present the dromion type localized 
solutions of eq.(l), the so-called induced localized structures/or induced dromions[14] for the 
potential u(x,y,t). This is possible by utilising the freedom in the choice of the arbitrary 
functions m lfi and m\j. For example, if we make the ansatz 

mu(p) = Km 1R (p) = tanh(p R ), (34) 

u = 2t7(£ — rjKjsech 2 p R sech[r]x + tanhp R — r)x ], (35) 

where p R = y — £i and A; is a real constant. Similarly, the expressions for the spin can 
be obtained from eqs.(33). The solution (35) for u{x,y,t) decays exponentially in all the 

— * 

directions, eventhough the spin S itself is not fully localized though bounded. Analogously 
we can construct another type of "induced dromion" solution with the choice 

dpR 

m u = Km 1R 



(PR + Po) 2 + 1 

where p an d m are constants, so that 



= / 7 , P \ 2 , ; + mo, (36) 

J (0R + 0nV + l 



2r](^ — ny) , 2 
u(x,y,t) = ■ r^— -sech 

[PR + PO) 2 + 1 



, f d PR 

v x + — — , : N9 , -i - v x o 



(37) 



(PR + Po) 2 + 1 

Generalizations of these solutions are also possible, which will be considered elsewhere. 

Finally, we note here that we have a non-isospectral problem, as the spectral parameter 
A satisfies eq.(15). The above presented solutions correspond to the constant solution of 
eq.(15), that is A = A x = constant. One may consider other interesting solutions of eq.(15). 
For example, one can have a special solution 
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A = Ai = rj(y, t) + i£(y, t) = y + k _+^ : (38) 

where b, k and r] are real constants. Corresponding to this case, we may call the solutions of 
eqs. (1) and (12) as breaking solitons[19]. Using the Hirota method, one can also construct 
the breaking 1-SS of eq.(l) associated with (38). For this purpose, we take gi in the form 

g = 9l = eX p x, x = ax + m + c = xr + iXi, (39) 

where a = a(y,t), m = m(y,t) and c = c(t) are functions to be determined. Substituting 
(39) into the first of eq.(28a), we get 

ia t + aa y = 0, im t + am y = 0, iA t + Aa y = 0, (40) 
where A = exp(c). Particular solutions of eqs. (40) have the forms 

r]-i(y + k) (y + k + irA A 

where 77, k, b and A are some constants. From eqs. (28 b,c), we obtain 

B jM lX-T r - (42) 

Now, we can write the breaking 1-SS of eq.(l) (using equations (25), (39)-(42)) 

. 2r]expi(xi + (j))(y + k-ir])[(y + k)coshz~ir]sinhz} 

S^{x,y,t) = 3 — , (43a) 

[(y + k) 2 + ri 2 }* cosh?z 

S 3 (x, y,t) = l- ^ y + k)2 + r]2] sech 2 z, (436) 

where z = jfex - \ In [(y + kf - V 2 ] + ^ ^ = In | |, X i = -^x + 

m i (^zr 2 )- We see that the solution (43) corresponds to an algebraically decaying solution 
for large x, y. 

Finally, we note that eq.(l) is a particular case of the following family of (2+1) dimen- 
sional equation 
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S t = {S AS y + uS} x + F, (44) 

where u and F satisfy the eq.(lb) and S-F = respectively. Eq.(44) admits many integrable 
reductions, for example, 

a) the isotropic M-I eq.(l), when F = 0; 

— # — * — * 

b) the anisotropic M-I equation, when F = S A AS, where A = diag{a\, a<i, CI3), and 

c) the M-II equation, when F = mvS x + nS y , where v x = k(S%) y , m, k, n are constants; 
and so on. All of these equations are integrable in the sense that the corresponding Lax 
representations exist [7] and their gauge equivalent counterparts can be constructed[10](see 
also [20]- [21]). So further studies of them will give more insight into the structure of nonlinear 
spin excitations in (2+1) dimensions. 
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